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Recently the Fermi-LAT data have revealed two gamma-ray emitting bubble-shaped structures 
at the Galactic center. If the observed gamma rays have hadronic origin (collisions of accelerated 
protons), the bubbles must emit high energy neutrinos as well. This new, Galactic, neutrino flux 
should trace the gamma ray emission in spectrum and spatial extent. Its highest energy part, above 
20-50 TeV, is observable at a kilometer scale detector in the northern hemisphere, such as the 
planned KM3NeT, while interesting constraints on it could be obtained by the IceCube Neutrino 
Observatory at the South pole. The detection or exclusion of neutrinos from the Fermi bubbles will 
discriminate between hadronic and leptonic models, thus bringing unique information on the still 
mysterious origin of these objects and on the time scale of their formation. 

PACS numbers: 95.85.Ry, 14.60.Pq, 98.70.Sa 



Gamma-ray data in the 1-100 GeV range from the 
Fermi-LAT pQ show a new, unexpected, feature of our 
galaxy: two huge spheroidal structures (Fermi "bubbles" , 
FB from here on), extending up to ~ 8 — 9 kpc (50°) out 
of the galactic center on either side of the galactic disk [2j 
[3]. Intriguingly, the bubbles coincide spatially with the 
WMAP "haze" in microwave [4] and the thermal X-ray 
emission seen by ROSAT [5|. The origins of these three 
phenomena, and whether they are related, are matters of 
intense debate 0E1EH9]. 

If we limit ourselves to the Standard Model of parti- 
cle physics, there are two main models for the FB. One 
is the leptonic origin: Compton scattering of relativistic 
electrons on photons (microwave, or optical/ultraviolet), 
where the electrons originate from shocks in the outflow 
of the Galactic center black hole, Sgr A* [3 , 8 , or from 
an episodic activity of Sgr A* due to the capture of a star 
[9] or to a localized star formation (SF) event [10 j. The 
second scenario is the hadronic origin: collision of accel- 
erated protons on background protons in the bubble gas 
[7]. The collisions produce 7r° mesons, which then decay 
as 7T° — ^ 77. Here the accelerated protons originate in 
the supernova remnants (SNRs) near the Galactic center. 
The two scenarios differ widely in the cooling times of the 
accelerated particles, and therefore in the predicted age 
of the bubbles: millions of years for the leptonic model, 
versus billions for the hadronic case. Thus, the origin 
of the FB will teach us about the time scale of Sgr A*'s 
activity and about the Galactic star-formation history. 

A major signature of the hadronic models - and there- 
fore a discriminator between the leptonic and hadronic 
models of the FB [7J [11] - is the presence of a high 
energy, ~ 1-10 TeV neutrino (v) counterpart of the 
7 rays (see, e.g., [E]). Indeed, pp collisions produce 
charged pions, 7r ± , as efficiently as 7r°, and muon and 
electron neutrinos (z/^,z/ e ) are then produced by decay, 
e.g., 7r + — >• ji +J r y^^y^i +e + + v e -\- v^. This mechanism 



is also the origin of the observed flux of neutrinos from 
cosmic-ray interactions in the Earth's atmosphere. 

In this Letter we present the first study of neutrinos 
from the FB, in the hadronic model, and show that they 
are detectable at a kilometer scale experiment, like the 
existing IceCube [13 , or the planned KM3NeT [14], de- 
pending on the detector's location. This is a new signal, 
whose detection or exclusion will contribute to under- 
standing the physics of the bubbles. 
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FIG. 1: The contour of the Fermi bubbles in equatorial 
coordinates (shaded), adapted from 3 . Considering Z as 
the zenith angle with respect to a neutrino telescope at a 
latitude +43° (KM3NeT), the area below each curve has 
cos Z < —0.66 at a given time of the day (labels on curves). 

Figure [I] shows a sketch of the FB in equatorial co- 
ordinates. They face the Earth's southern hemisphere, 
extending to ~ 60° below the equator, and subtend a 
solid angle ^5 = 0.808 sr [3]. The bubbles are uniformly 
bright in 7 rays, with sharp edges and a hard spectrum, 
dN 1 /dE oc E~ 2 [3]. 

In the hadronic model of Crocker and Aharonian 
(CA) [7 , the bubbles are powered by SNRs, which are 
widely believed to be the factories of cosmic rays up 
to 10 15 -10 17 eV [17 . The low-density bubble interiors 
(average density (uh) ~ 10 -2 cm -3 [3]) are created by 
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prolonged SF activity near the Galactic center, which 
forms a high-velocity bipolar wind. This wind carries 
the cosmic ray protons produced in the SNRs to fill the 
bubble cavities. The interactions of these protons with 
the diluted hot gas produce the FB as observed by the 
Fermi-h AT . The associated time scale for proton cooling 
is ~ 5 x 10 9 yr [7], and therefore a SF period of at least 
this duration is required in the CA model. 

To calculate the 7 and v fluxes from pp interactions in 
the FB we use an injected proton spectrum in the bubble 
cavities, at steady-state, of the form 



N p (E) = N E- k eM-E/E ), 



(1) 



where TVo, k and £"0 are the normalization factor, spectral 
index and cutoff energy respectively. The cutoff energy, 
Eq °° 1 — 10 PeV, is motivated by cosmic-ray accelera- 
tion in SNRs [17 , and k ~ 2 is expected in a Fermi ac- 
celeration mechanism. The index of the injected protons 
remains hard, in the FB, due to a saturation condition 
realized when the acceleration (in the SNRs) time, pp en- 
ergy loss (in the FB) time and escape (from the FB) time 
satisfy a relation: t acc < ti oss < t esc [7 . Using a param- 
eterization from the Monte Carlo code SIBYLL [T8] , we 
obtain the fluxes of gamma rays and neutrinos from 7r°, 
7r ± decay, above > 10 GeV. At lower energies, the fluxes 
can be modeled with a delta-function approximation. For 
example, the 7 flux is given by 



$7(^7) = 



4ttD 2 K, 



- r dE„ ^A^ N P (E C ). (2) 



Here D is the distance to the bubbles, K n « 0.17 is the 
mean fraction of proton kinetic energy converting to 7r°, 
Ejr^h — E 1 + ra^/4E 7 is the threshold pion energy, <j pp 
is the inelastic pp cross section, and E c = E^/K^ + m p . 
The parameter tp ~ 0(1), and is adjusted to match the 
E > 10 GeV Monte Carlo results. Note that there is a 
~ 10% uncertainty in the hadronic models [18 . We also 
neglected a contribution of ~ 10% or less (depending on 
the unknown magnetic field in the FB) due to Compt on- 
scattered photons off secondary e ± . Our calculation gives 
the fluxes of 7 and of v e , v e , in the expected flavor 

ratio \v e — 2 : 1 for both neutrinos and antineutrinos. 
Here we focus on for which detectors are the most 

sensitive. The fluxes of 7 and v ^ + v ^ are shown in 
Fig. [2j for k = 2.1. They are compared to the flux of 
atmospheric v^, which constitutes the main background. 
For this, we used the model in [19], extrapolated to fit 
the high energy IceCube data [16]. At our energies of 
interest, flavor oscillations are negligible for atmospheric 
neutrinos, due to the short propagation length. However, 
they are effective for the FB flux, so that the flux at 
Earth is about 1/2 of the unoscillated one (see e.g. [20 J. 

Let us comment on the energetics of the CA model. 
The total number of 7-ray producing cosmic-ray pro- 
tons, in steady state, in the FB is J dE N P (E) w 
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FIG. 2: The fluxes of 7 and of (total of neutrinos and 
antineutrinos) from the FB, for different proton cutoff energy 
Eo (Eq. |l])), compared with the atmospheric flux and the 
Fermi-LAT data 3 . Errors are la. 



4.3 x 10 57 ((n#)/10- 2 cm- 3 )- 1 ^^ kpc) 2 , where 
D — 8.5 kpc is the distance to the bubble cen- 
ters [3 j. The total energy in steady state cosmic- 
ray proton population is J dE N p (E)E w 5.4 x 
10 55 ((n#)/10- 2 cm- 3 )- 1 ^/^ kpc) 2 ergs. For com- 
parison, with a 2 x 10 37 erg s _1 power in 7 rays from 
both bubbles [26] and a ~ 5 x 10 9 yr lifetime of the FB, 
the total energy emitted in 7 rays is ~ 3 x 10 54 ergs. 
This is ~ 6% of the total cosmic-ray proton energy, an 
efficiency typical of hadronic models. 

In order to produce the observed uniform surface 
brightness of FB, a volume emissivity of the form (R 2 — 
r 2 ) -1 / 2 ; (r < R ~ 3.5 kpc) seems required [3]. Whether 
this can in principle be obtained in the hadronic scenario 
has not been studied and is beyond our purpose here. A 
possibility is that anisotropics can arise from the diffu- 
sion of cosmic rays in the magnetic field generated by the 
shocked shells of the bubble [27] . In any case, the v flux 
should trace the 7-ray flux from the FB, and therefore 
have the same uniform projected intensity. 

Under this assumption, we can write the number of 
events accumulated in a certain exposure time t expi and 
from a certain zenith bin [#i,02]> (with respect to the 
detector's position) as follows: 

pt exp p poo 

N v = dt dn / dE $(E)A U (E,6 Z ) 

Jo J ei<e z <e 2 jEth 

~t exp (n(6 1 ,6 2 )) t dE *(E){A v (E)) e , (3) 

^ Eth 



with a completely analogous expression for . Here A v is 
the detector's effective area, Z is the zenith angle with re- 
spect to the detector, and £(£) indicates that the integral 
in the solid angle is done over the region of the bubbles 
for which the condition on the zenith angle is satisfied. 
This region depends on the time of the day, due to the 
revolution of the Earth (Fig. [TJ. The approximation in 
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Eq. ([3| uses the time- averaged solid angle, (£}(0i,02))t? 
and the effective area averaged over the zenith bin of in- 
terest, (A u (E))q. The latter is adequate for sufficiently 
narrow bins, since A v varies gently with Z . 

As a benchmark, we use twice the effective area of the 
IC40 configuration of IceCube [16!, A v = 2Aic4o, which 
well represents the full IceCube configuration [21] and is 
reasonable for KM3NeT. We take a 10-year exposure, and 
consider only the v ^ + v ^ flux in sub- horizon directions, 
6 Z > 90°, as these channels dominate the effective area. 

Under this condition, a detector in the northern hemi- 
sphere is more sensitive to the FB compared to a south- 
ern one since the Galactic center is in the southern hemi- 
sphere (Fig. [TJ. This appears clearly in Fig. [3j which 
shows the ratio 02))t/£tb, f° r different zenith bins 

of interest, and as a function of the detector's latitude, 
0i at . IceCube, at Q\ at = —90°, is sensitive only to the 
small portion of the upper bubble that lies in the north- 
ern hemisphere (Fig. [TJ, for which (Q) t « 0.02 sr. For 
KM3NeT we used lat = +43°, the position of the cur- 
rent ANTARES detector in the Mediterranean [22 , and 
we have a much more encouraging total of (Q(cosO z < 
0))t ~ 0.57 sr (Fig. [3j. Note that this location is nearly 
optimal (maximum solid angle) for events in the deepest 
bin, cos0 z < —0.66, where the signal to background ratio 
is the highest (see below). This bin corresponds to a cone 
of opening angle ~ 48° , which is wide enough to contain 
both bubbles almost completely when the center of the 
lower bubble is at Z ~ 180° from the detector (Fig. [T]). 
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FIG. 3: Daily- average fraction of the solid angle subtended by 
the Fermi bubbles, = 0.808 sr, that falls in a given bin of 
zenith angle, Z (with respect to the detector), as a function 
of the detector latitude. The bins in cos0 z are: [—1,-0.66] 
(long dashed), [-0.66,-0.33] (short dashed), [-0.33,0] (dot- 
ted) and [-1,0] (total sub-horizon, solid). 

Figure [4] (top panel) shows the number of signal and 
background events at KM3NeT, as a function of E t h for 
the entire region of sensitivity (cos# z < 0) and for the 
deepest zenith bin (cos# z < —0.66). They are in rough 
agreement with the earlier estimates in [7]. 

Counting on the good directional sensitivity of the 
detector, we restrict the analysis of the background to 



the daily- averaged solid angles subtended by the bub- 
bles, for better signal discrimination. For the entire 
sub- horizon solid angle (solid line in Fig. |3|, and for 
Eq = 3, 10 PeV, the signal exceeds the background above 
~ 50 — 100 TeV, where ~ 20 — 100 signal events are ex- 
pected. We found that Eth ~ 20 — 30 TeV gives the 
maximum statistical excess of the signal over the back- 
ground, A = N 3 ig/ ^/N S ig + N^g, which is between 4 
and 6 a (i.e., A ~ 4 — 6). 

For the deepest zenith bin the event rate falls faster 
with Eth , compared to the entire region of zenith, reflect- 
ing the stronger opacity of the Earth to z/s crossing its di- 
ameter. However the atmospheric flux is also suppressed 
due to the strongest absorption of muons. Therefore the 
signal is higher than background already above ~ 20 TeV 
in the most optimistic case. The statistical significance 
of the signal is lower due to the lower number of events, 
but it exceeds 3a for Eq = 3, 10 PeV, nevertheless. 

The zenith distributions of the signal and background 
are given i n Fig, [j] (bo ttom panel), together with the 
quantity ^N sig + N bkg in the form of error bars. In all 
cases, the numbers of signal events in the three bins are 
similar. This is because the strongest neutrino absorption 
in the deepest bin (lower A u ) compensates for the larger 
solid angle. The figure also confirms how the number of 
events in the deepest bin and in the total sub-horizon 
region are the best indicators of the effect, in terms of 
statistical significance, while the other two bins individ- 
ually suffer from low statistics and high background. 

Some remarks are in order. We have taken k = 2.1, but 
results obtained with a much more conservative cosmic- 
ray proton spectrum, k = 2.3, also fit the 7-ray data [7]. 
In this case, the v event rate is dominated by background, 
and A ~ 1.2^/t exp /10 yr, for cos0 z < 0- A ~ 6 times 
longer exposure would be required to reach 3a. Note, 
however, that using A > 3 as criterion of statistical sig- 
nificance is conservative. Indeed, an excess below 3a in 
number of events is likely to be further substantiated, in a 
full fit of the data, by other indicators like the correlation 
with the position and extent of the FB. 

Briefly, let us comment on the sensitivity of a south 
pole detector. Due to the suppression in the solid angle 
and the near-horizon direction of the signal, the num- 
ber of signal events is low and the background level is 
enhanced compared to KM3NeT. For the best case sce- 
nario of E = 10 PeV, we find N sig = 8.0 and N bkg = 59 
for E th = 10 TeV, and N sig = 2.4 and N bkg = 1.7 for 
E t h = 100 TeV. Thus, an excess above background will 
not be statistically significant, unless an exposure at least 
~ 6-9 times larger is achieved or efforts to gain sensitivity 
to events above the horizon [23] and to better discrim- 
inate background (e.g., by detecting tau neutrinos [24] ) 
are successful. While a discovery with IceCube seems un- 
likely, still an analysis of the IceCube data might place 
interesting constraints on the parameters: for example, 
it may be able to exclude larger, less natural, values of 
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FIG. 4: Number of signal and background events in 10 years, 
at a km scale detector at +43° latitude, for different cutoff 
energy Eq in the proton spectrum. Top: as a function of 
the energy threshold, E t h- Bottom: zenith distribution, for 
E t h = 25 TeV. Expected la error bars are shown. 



Eo or harder proton spectra, k < 2.1. 

In summary, the FB are a potential new source of 
galactic high energy neutrinos. For spectral cutoff Eo ~ 
1-10 PeV, and spectral index k ~ 2.1 for cosmic-ray pro- 
tons - motivated by particle acceleration in SNRs - the 
neutrino flux could exceed the atmospheric background 
above ~ 20 — 50 TeV and be seen at km-scale detectors. 
For a steeper spectrum, detection will require longer ex- 
posures or larger detectors. IceCube is strongly disad- 
vantaged by its southern location, but it could obtain 
important constraints, and might achieve detection with 
an increased sensitivity to the sky above the horizon. Our 
results motivate efforts in this direction. For a northern 
detector, observing the FB with a decade of operation is 
a realistic goal, that contributes to make the case for the 
planned KM3NeT in the Mediterranean. It is encourag- 
ing that the Italian node of KM3NeT is funded [25]. 

The observation of neutrinos from the FB would 
strongly support the hadronic origin of 7 rays from these 
mysterious objects and favor the model for ~ 10 9 years 
scale activity at the Galactic center to form the FB 0. 
The shape of the neutrino spectrum would help to dis- 
tinguish this scenario from more exotic neutrino-emitting 
explanations of the FB, like dark matter annihilation or 



decay [6 j. A null result would indicate a softer proton 
spectra or a non-hadronic mechanism as the primary ori- 
gin of 7 rays and favor ~ 10 6 years scale activity by Sgr 
A* to form the FB [3]. 

As neutrino and 7-ray detectors progress, it is likely 
that their interplay will ultimately identify the nature 
and origin of the FB. Implications will be far reaching, 
showing what high energy phenomena are at play in our 
galaxy and, by extension, in other galaxies too. 
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Fermi Bubble y-ray data 
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FIG. S-l: Contribution to the observed 7-ray spectrum from 
Compton-scattered CMB photons, labeled SIBYLL (IC only), 
by ^ decay in case the bubble magnetic field is 3 fiG. 
The IC component, the normalization of which depends on 
the unknown bubble magnetic field, is ignored in Fig. [2] of 
the main text. Also shown is the 7-ray spectrum from an- 
other hadronic model, labeled QGSJET (tv° only), for the 
same injected proton spectrum, to illustrate hadronic model 
uncertainties. 
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FIG. S-3: Number of signal and background events in 10 
years, for 3 zenith- angle bins and the total, in a km scale 
detector at +43 degree latitude (KM3NeT). The events are 
calculated as a function of the neutrino threshold energy and 
for up to 1 PeV energy. The event rates in the top left and 
bottom right panels are reported in Fig. ^] of the main text. 
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FIG. S-2: Neutrino and antineutrino flux components before 
oscillation from ^ chain decay for the total + flux, 
labeled as Total zy^, plotted here and in Fig. [2] of the main 
text. The z^ e +z> e component is very similar to the tv — » ji —> 
component plotted here. 



FIG. S-4: Statistical significance of the bubble neutrino signal 
for the events in Fig. |S-3| above. 
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TABLE S-I: Number of signal and background events for a 
km scale detector at the South pole (IceCube). The daily- 
averaged solid angle is 0.02 sr (see Fig. [3] in the main text). 
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FIG. S-5: Gamma-ray and neutrino flux models from the 
Fermi bubbles in case of a conservative cosmic-ray proton 
spectral index of —2.3, instead of —2.1 used in Fig. [2] of the 
main text and all previous analyses. 
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FIG. S-6: Detection prospects of neutrinos from the Fermi 
bubb les f or the cosmic-ray proton spectral index of —2.3 
(Fig. |S-5[ ), at a km scale detector at +43 degree latitude 
(KM3NeT). Top - Number of signal and background events 
in 10 years. Bottom - Detection significance for 10 years 



